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^ ' Abstract 

^' We investigate the mass spectrum of Nucleoli and Delta (and its 

^I^ • counterparts with strange and charm), and their excited states, in quenched 

lattice QCD with exact chiral symmetry. For each light baryon, we use 
K^ I 23 masses to determine the coefficients of the mass formula in quenched 

^ ■ chiral perturbation theory. By chiral extrapolation to m.,^ = 135 MeV, 

<^ ■ we obtain Mjv = 958(26) MeV, Mn* = 1553(42) MeV, Ma = 1216(32) 

MeV and Ma* = 1611(17) MeV, which are identified with A^(939)Pii, 
iV(1535)S'ii, A(1232)P33 and A(1620)S'3i respectively Further, we di- 
rectly measure the masses of fi", Mq = 1648(60) MeV, and its excited 
state, Mfi* = 1935(48) MeV; as well as the triply charmed baryon r2++, 
M^++ = 4931(22) MeV, and its excited state, M^++* = 5185(35) MeV. 
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One of the objectives of lattice QCD is to compute the hadron masses 
nonperturbatively from the first principles. For hadrons composed of charm 
and strange quarks (i.e., without u, d light quarks), their masses can be directly 
measured on presently accessible lattices. However, for hadrons containing m, d 
light quarks, the performance of the present generation of computers is still 
quite remote from what is required for computing their masses at the physical 
scale (e.g., m.,^ ~ 135 MeV), on a lattice with enough sites in each direction 
such that the discretization errors as well as the finite volume effects are both 
negligible comparing to the statistical ones. Nevertheless, even with lattices of 
moderate sizes, lattice QCD can determine the values of the parameters in the 
hadron mass formulas of the (quenched) chiral perturbation theory. Then one 
can use these formulas to evaluate the hadron masses at the physical scale, as 
well as to determine the quark masses. 

In this paper, we study spin 1/2 and spin 3/2 baryons containing three 
quarks of the same mass^, and obtain their masses for 30 quark masses in the 
range 0.03 < iriqa < 0.8 (i.e., from ms/2 to r/ic)- 

For spin 3/2 baryons composed of three strange quarks (f^^), and three 
charm quarks (fi^i^^), as well as their excited states, we can measure their 
masses directly on the lattice. Here the strange quark bare mass {rrisa = 0.06) 
and charm quark bare mass {rrica = 0.8) are fixed by requiring the masses 
of the corresponding vector mesons to agree with the masses of 0(1020) and 
J/ V- (3097). 

For baryons composed of light quarks, we use 23 masses to determine the 
coefficients in the baryon mass formula in quenched chiral perturbation theory 

HI 

M = Co + cim^ + ciml In(m^) + C2m^ + c^ml + 0{ml In m^) (1) 

Then, extrapolating to the physical pion mass m^ = 135 MeV, we obtain the 
masses of A^(l/2+), A^*(l/2-), A(3/2+), and A*(l/2-) respectively. For our 
data, we find that we have to set ci = 0, otherwise the errors of some coef- 
ficients would become unacceptable (> 100%). For A*, we also have to set 
C3 = 0. Evidently, if one can measure more baryon masses in the chiral regime 
(?7i^ ^ 2v2vr/^), then the coefficients can be determined more precisely. How- 
ever, one should limit mT^L > 3 ~ 4, otherwise the baryon masses would suffer 
from finite size effects, which in turn would lower the reliability of the resulting 
coefficients. For our lattice of size 20^ x 40 at /3 = 6.1, we find that ruqa > 0.03 
gives rrLj^L > 4 in the spatial direction, which should be sufficient to suppress 
the finite size effects. 

The interpolating operators for Nucleon and A can be written as 

Nxa = ^abc^xabi^'y^) a/sd-x/Sc'^xaa (2) 



"'^In this paper, we work in the isospin hniit, rn^ — m^. 
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and 

^xa■,^l = eabcd-xabi^l5lfj.)a/3d-xf3cd-xaa , /i = 1, 2, 3 (3) 

where u and d denote the quark fields, the superscript T denotes the transpose 
of the Dirac spinor, C is the charge conjugation operator, eatc is the completely 
antisymmetric tensor, and x, {a, b, c} and {a, j3} denote the lattice site, color, 
and Dirac indices respectively. For A(A*), we average over fi = 1,2,3 to 
increase the statistics. Note that the "diquark" operator in Q transforms like 
a vector, thus the even parity state of Q can overlap with A (3/2"*"), in which 
the quarks are in S-wave; while its odd parity state can overlap with A*(l/2~), 
in which the quarks are in P-wave. 

Now it is straightforward to work out the baryon propagator in terms of 
quark propagators. In lattice QCD with exact chiral symmetry, quark propa- 
gator with bare mass rrig is of the form [Dc + mq)~^ [3], where Dc is exactly 
chirally symmetric at finite lattice spacing. In the continuum limit, {Dc+mq)~^ 
reproduces [jfj.id^ + ^^4^) + mq\~^. For optimal domain- wall fermion [2] with 
Ng + 2 sites in the fifth dimension, 

Dc = 2mo: 



S{H, 



i-nf=iT, 



- s 



1 — LUgHw 
1 + LOsH^,, ■ 



Hw = l^D^ 



where Dyj is the standard Wilson Dirac operator plus a negative parameter 
—mo (0 < rriQ < 2), and {us} are a set of weights specified by an exact formula 
such that Dc possesses the optimal chiral symmetry |2]. Since 

{Dc + m^)"^ = (1 - rmq)~^[D^^{mq) - r], r 



2r7io 



where 

D{'mq) = mg + (mo - mg/2)[l + 'y^S^H^)], 

thus the quark propagator can be obtained by solving the system D{mq)Y = 1 
with nested conjugate gradient |1|, which turns out to be highly efficient (in 
terms of the precision of chirality versus CPU time and memory storage) if 
the inner conjugate gradient loop is iterated with Neuberger's double pass 
algorithm [3]; . For more details of our scheme of computing quark propagators, 
see Ref. j6]. 

We generate 100 gauge configurations with Wilson gauge action at /5 = 6.1 
on the 20^x40 lattice. Fixing m-o = 1.3, we project out 16 low-lying eigenmodes 



of \H.w\ and perform the nested conjugate gradient in the complement of the 
vector space spanned by these eigenmodes. For A^^ = 128, the weights {us} 
are fixed with X^in = 0.18 and Xmax = 6.3, where Xmin < H\Hw\) < >^max 
for all gauge configurations. For each configuration, quark propagators are 
computed for 30 bare quark masses in the range 0.03 < rriqa < 0.8, with 
stopping criteria 10~^^ and 2 x 10^^^ for the outer and inner conjugate gradient 
loops respectively. Then the norm of the residual vector of each column of the 
quark propagator is less than 2 x 10"^^ 



-11 
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\\(D^ + mg)Y - I|| < 2 X 10"^ 
and the chiral symmetry breaking due to finite A^^ is less than 10~^^, 



yty 



-1 



< 10 



for every iteration of the nested conjugate gradient. 

After the quark propagators have been computed, we first measure the 
pion propagator and its time correlation function, and extract the pion mass 
{mT,a) and the pion decay constant {fnd). With the experimental input /^ = 
132 MeV, we determine the inverse lattice spacing a^^ = 2.21(3) GeV. (Our 
procedure has been described in Ref. E].) 

Next we compute the nucleon propagator {N^aNyi) and its time correlation 
function C{t). Then the average of C{t) over gauge configurations is fitted by 
the usual formula 



i^^(Z+e-'"+'^* - z_e-"^-"(^-*)) + L^(^^e-™+"(^-*) - Z.e 



m-at\ 



where m± are the masses of even and odd parity states. Thus, one can use 
parity projector (1 ± 74) /2 to project out two amplitudes. 

Now the problem is how to extract m± from A± respectively. Obviously, for 
sufficiently large T, there exists a range of t such that, in A±, the contributions 
due to the opposite parity state are negligible. Then m^ can be extracted by 
a single exponential fit to A±, for the range of t in which the effective mass 
rrieffit) = ln{A-i-{t) / A±{t + 1)) attains a plateau. On the other hand, if T is 
not so large, then it may turn out that the heavier mass, say m_ (assuming 
?Ti_ > «!+), could not be easily extracted from A_ due to the non- negligible 
contributions of the (lowest lying) even parity state. Further, if T is too small, 
then one even has difficulties to extract the mass of the lowest lying state. For 
our lattice with T = 40, it is sufficiently large to extract both m+ and m_ 
from y4_|_ and A_ respectively. 
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Figure 1: The nucleoli masses versus the pion mass square, for N{J^ = 1/2"'') 
and N*{J^ = l/2~) respectively. The solid lines are fits to (Q) with Ci = 0. 



In Fig. n the nucleon masses of even and odd parity states are plotted 
versus {mT^aY, for m^r smaller than the chiral cutoff A^ = 2v^7r/^. For the 
J^ = 1/2+ ground state nucleon, the 23 nucleon masses can be fitted by 

Mnu = 0.417(13) - 1.350(168) x (m^a)2ln(m^a) + 0.727(134) x {m^af 

+0.458(144) X {m^af (4) 

with x^/d.o.f. < 0.003. At the physical pion mass 171^^ = 135 MeV (the y-axis), 
dH) gives Mn = 958(26) MeV, which is naturally identified with A^(939)Pii. 
For the J^ = l/2~ excited nucleon, the 23 masses can be fitted by 

MN*a = 0.689(21) -0.949(291) x (m^a)Mn(m^a) + 0.904(236) x {m^af 
+0.211(250) X {m^af (5) 

with xVd.o.f. < 0.007. At m^ = 135 MeV, © gives M^* = 1553(42) MeV, 
which is identified with A^(1535)5'ii. 




Figure 2: The nucleon mass square versus the pion mass square, for N{J^ 
1/2+) and N*{J^ = 1/2") respectively. 



At this point, it is instructive to examine whether our data of nucleon mass 
can accomodate the linear term citti^. To this end, we take the square of the 
nucleon mass, then the higher order terms are further suppressed, and the 
linear term (if any) would become more prominent, especially at small rriTrCi- 



In Fig. we plot (Matc)^ and {Misf*a)'^ versus 
well fitted by 



m^ra 



P, which turn out to be 



(Mtvo)^ = 0.196(9) + 2.388(55) x {m^aY 
(MN-af = 0.5174(192) + 2.829(124) x {m^af 



(6) 
(7) 



Obviously, the term CiJ(m.,^a)'^ does not exist in our data of nucleon masses, 
otherwise, we would see deviations from the fitted lines at small {mT^a). 

Note that © and © give Mn = 1001(22) MeV and M*^ = 1607(30) 
MeV at ttItt = 135 MeV, in good agreement with the masses of iV(939)Pii 
and A^(1535)S'ii. This seems to justify the ansatz, namely, if the number of 



baryon mass data points is too small (e.g., less than 10) and/or the baryon 
masses are measured only at large quark masses, then they may not feasible for 
fitting the mass formula ([Q), however, they still could be used for the linear fit: 
{Ma^ = A + B{mT^a)'^, which may turn out to provide a good extrapolation 
to the physical rrin- 

For recent quenched lattice QCD studies of A^ and A^*, see Refs. [Z]-[I21- 
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Figure 3: The A masses versus the pion mass square, for A{J^ = 3/2+) and 
A*(J^ _ i/2~) respectively. The solid lines are fits to (jT]), with ci = for A, 
and with ci = C3 = for A*. 



Next we compute the time correlation function of A, and extract the masses 
of even and odd parity states. 

In Fig. |3[ we plot the masses of A versus m^ for even and odd parity states 
respectively. The lowest lying state has even parity, which implies that it has 
J^ = 3/2"*" since its three quarks are in S-wave, and symmetric in spin. Using 
dD) with ci = to fit the 23 masses of A, we obtain 
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Figure 4: The Delta mass square versus the pion mass square, for A(J^ 
3/2"*") and A*{J^ = l/2~) respectively. The solid lines are linear fits. 



M^a = 0.530(16) - 1.685(218) x (m^a)Mn(m^a) + 0.552(176) x {m^af 

+0.682(187) X {m^af (j" = 3/2+) (8) 

At m^ = 135 MeV, © gives Ma = 1216(32) MeV, which is identified with 

A(1232)P33. 

For the negative parity state A*, its quarks are in P-wave, thus its J^ = 
1/2^. Using (^ with ci = C3 = to fit the 23 masses of A*, we obtain 

MA-a = 0.717(8) - 0.724(34) x (m^aY ln{m^a) + 1.156(12) x {m^aY (9) 

At m^ = 135 MeV, © gives Ma* = 1611(17) MeV, which is identified 
with A(1620)5'3i. This is the first lattice QCD determination of the mass 
of A(1620)53i. 



Again, we check the ansatz of baryon mass square versus pion mass square. 
In Fig. m we plot (Mao)^ and (MA*a)^ versus {mT,aY. They can be fitted by 
the straight fines: 

(M^af = 0.330(13) + 3.078(87) x {m^af [j" = 3/2+] (10) 

(MA*a)2 = 0.536(17) + 3.000(109) x {m^ay [J^ = 1/2"] (11) 

At m^ = 135 MeV, tfiey give Ma = 1292(25) MeV and M^, = 1635(26) MeV, 
in good agreement A(1232)P33 and A(1620)S'3i respectively. 

Finally, we turn to tfie fieavy baryons wfiicfi can be measured directly on 
our lattice. Now replacing tfie down quark in A~ witfi tfie strange quark, we 
fiave tfie ^2". Similarly, replacing tfie up quark in A"*""'" witfi tfie cfiarm quark, 
we obtain tfie triply cfiarmed baryon f2++(J^ = 3/2+), wfiicfi so far fias not 
been discovered in fiigfi energy experiments. 

To determine tfie bare masses of strange quark and cfiarm quark, we extract 
tfie mass of vector meson from tfie time correlation function 

Cp{t) = ^ E E tr{7/.(^c + m,)-i(0, xh,{D, + m,)-\x, 0)} 

At ruqa = 0.06, tfie vector meson fias mass Ma = 0.4638(32), wfiicfi gives 
M = 1025(7) MeV, in good agreement witfi tfie mass of 0(1020). Tfius, at 
rriqa = 0.06, A~ becomes Q~ , and its mass can be extracted from tfie time 
correlation function. Our result is MQa = 0.746(27), wfiicfi yields 

Mn- = 1648(60)MeV (12) 

in good agreement witfi tfie mass of Q^{1Q72). 

Similarly, at rriqa = 0.80, tfie vector meson fias mass Ma = 1.3830(15), 
wfiicfi gives M = 3056(3) MeV, in good agreement witfi tfie mass of J/\E'(3097). 
Tfius, at rUqa = 0.80, A++ actually is f2^+, tfie triply cfiarmed baryon witfi 
J^ = 3/2+, wfiicfi fias not been observed in experiments. However, we can 
determine its mass from tfie time correlation function, and our result is Ma = 
2.23(1), wfiicfi gives 

Mf^++ = 493 1(22) MeV [quark content=(ccc), j" = 3/2+] (13) 

Also, for tfie corresponding J^ = l/2~ excited baryons, our results are 

Mo* = 1935(48)MeV [quark content=(sss), j" = 1/2^] (14) 

M^++. = 5185(35)MeV [quark content=(ccc), j" = 1/2"] (15) 

Note tfiat (fT^ - (fT3j) are tfie first lattice QCD predictions for tfie masses of tfiese 
baryons. Even tfiougfi we fiave not estimated tfie errors due to quencfied ap- 
proximation, discretization, and finite lattice size, we suspect tfiat tfie resulting 
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error in each mass of (fT^ - (fT3jl is less than 80 MeV. A possible candidate for 
the n* in (HH) could be the ^"(2250) ^. 

Finally we transcribe the bare masses of strange quark {mga = 0.06) and 
charm quark [rrica = 0.80) to MS at /i = 2 GeV. Using the lattice renormal- 
ization constant Zm = Z^^ (where Zg is the renormalization constant for i/jip) 
to one loop order |14J, we obtain 



m 
m 



^^(2 GeV) = 114(9) MeV (16) 

J^^(2 GeV) = 1520(38) MeV (17) 



where the errors are estimates based on the errors of the vector meson masses 
with respect to 0(1020) and J/^(3097). 
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David Lin for discussions. 
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